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ABSTRACT 

In this first paper we discuss the linear theory and the background evolution of a new 
class of models we dub SCDEW: Strongly Coupled DE, plus WDM. In these models, WDM 
dominates today’s matter density; like baryons, WDM is uncoupled. Dark Energy is a scalar 
field $; its coupling to ancillary CDM, whose today’s density is <C! 1 %, is an essential model 
feature. Such coupling, in fact, allows the formation of cosmic structures, in spite of very low 
WDM particle masses 100 eV). SCDEW models yields Cosmic Microwave Background 
and linear Large Scale features substantially undistinguishable from ACDM, but thanks to the 
very low WDM masses they strongly alleviate ACDM issues on small scales, as confirmed via 
numerical simulations in the II associated paper. Moreover SCDEW cosmologies significantly 
ease the coincidence and fine tuning problems of ACDM and, by using a field theory approach, 
we also outline possible links with inflationary models. We also discuss a possible fading 
of the coupling at low redshifts which prevents non linearities on the CDM component to 
cause computational problems. The (possible) low-z coupling suppression, its mechanism, 
and its consequences are however still open questions -not necessarily problems- for SCDEW 
models. The coupling intensity and the WDM particle mass, although being extra parameters 
in respect to ACDM , are found to be substantially constrained a priori so that, if SCDEW is 
the underlying cosmology, we expect most data to fit also ACDM predictions. 

Key words: large-scale structure of Universe, dark matter, dark energy - Galaxy: evolution, 
formation 


1 INTRODUCTION 

When SNIa Hubble diagrams revitalized ACDM, one could hardly 
guess how successful this cosmology would be in meeting cosmo¬ 
logical data. In spite of its coincidence and fine tuning problems, 
therefore, ACDM is surely the benchmark for any attempt to im¬ 
prove our understanding of cosmology. 

In this series of works we present a detailed analysis of a new 
class of cosmologies that preserves all successes of ACDM models 
on large and intermediate scales. In turn, they improve the agree¬ 
ment with_theob^rved_dark_jnatter_di^ributon on small scales 
(e.g. lOh et alj|201ll : I Walker & Penarrubiall201lh . at the same time 
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easing the coincidence and fine tuning problems of ACDM, while 
possibly launching a bridge between the reheating stages in the late 
inflationary regime and today’s Dark Energy (DE) nature. 

These cosmological models are based on the presence of three 
dark components: i) a very-warm dark matter component (WDM), 
uncoupled (like the baryons), and constituting the observable dark 
matter at the present time; ii) a scalar field 4> which acts as Dark 
Energy; iii) a peculiar Cold Dark Matter component (CDM), whose 
coupling to <0 played an essential role in cosmic history, but whose 
density at z = 0 naturally became almost negligible. We dubbed 
these models SCDEW (Strongly Coupled DE -b WDM). 


The basic features of these models were disc us sed i n 
two previous paper s (iBonometto. Sassi, & La Vaccal ( 12012h : 
iBonometto & Maininil (12014 ): BSLV12 and BM14 respectively, 
hereafter). BSLV12 dealt with background components: they 
showed that a purely kinetic scalar field $ strongly coupled to 
a CDM component would both exhibit densities oc a~* (a : 
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scale factor), during the radiative eras, expanding along a tracker 
solution, with (primeval) density parameters 0(0.1 %). The derel- 
ativisation of a further WDM component eventually turns these 
components into quintessential DE and a tiny CDM component. 
The radiative eras of these models are therefore characterized by 
three low density components, in top of the usual 7 ’s and u’s: 
WDM, coupled CDM and scalar field, accounting for constant 
fractions of the overall density, and sharing similar densities. 
BM14 then deals with fluctuation evolution, showing that CMB 
anisotropy and polarization spectra, in these models, are hardly 
distinguishable from the ACDM benchmark, while coupled CDM 
fluctuations continue to grow also between their entry in the 
horizon and matter-radiation equality, being so able, in spite of 
the low CDM density, to revitalize WDM fluctuations on scales 
suffering an early free streaming, even when the WDM particle 
mass rriw is quite small. 

The <1?-CDM coupling is therefore an essential feature in the 
early model evolution. After WDM has derelativized, instead, the 
<£>-CDM coupling is unessential, while leading to technical difficul¬ 
ties, as 5c (CDM fluctuation amplitudes) tend to approach unity, de¬ 
veloping then some early non-linearities, typically involving 1 % 
of the total mass. However, switching off the $-CDM coupling af¬ 
ter the break of “conformal invariance”, yields more comfortably 
tractable models. 

The plan of the paper is as follows. In the next Section we 
shall formulate the Lagrangian theory for strongly coupled Dark 
Energy cosmologies. It confirms the prediction on the coupled $- 
CDM component densities, while outlining possible links with in¬ 
flationary theories. In Section 3 we discuss the evolution of back¬ 
ground parameters, focusing on the exit from the primeval sta¬ 
tionary regime. We also widen the range of models in respect to 
BSLV12 and BM14 papers, by allowing for /) fading after the 
exit from conformally invariant expansion. In Section 4 we discuss 
fluctuation evolution, taking also into account a possible /3 fading. 
The effects of the early growth of CDM fluctuations are therefore 
analised, both as cause for the restart of WDM and baryon fluctu¬ 
ations, and for the possible formation of late CDM structures. We 
also show that, in order to mimic ACDM phenomenology through 
SCDEW models, the mass of WDM particles and the coupling are 
significantly constrained. In Section 5 fluctuation spectra are ob¬ 
tained and the model for N-body simulations of Paper II is con¬ 
veniently selected. The last Discussion Section outlines the predic¬ 
tions of SCDEW models, namely for what concerns their discrimi¬ 
nation from ACDM, whose main results are however faithfully re¬ 
produced. We also outline how SCDEW cosmologies, besides of 
easing low-scale ACDM conundrums, free us from its coincidence 
paradox and fine tunings. 


2 COUPLED DARK ENERGY IN THE EARLY 
UNIVERSE 

There are direct evidences that Dark Matter (DM) is a physical cos¬ 
mic component, mostly clustering with observable baryons. Doubts 
have been cast, on the contrary, on the true nature of Dark Energy 
(DE); here we assume it to be a self-interacting scalar field $, keep¬ 
ing essentially unclustered. In this scheme, both dark components 
interact with baryons and radiations just gravitationally, so that the 
energy pseudo-conseiwation equation 

= 0 (1) 


holds (here are the stress-energy tensors of DM and DE, 

their traces reading T'^'^’*)). This sum can vanish thanks to separate 
vanishings of both terms. The alternative option that 

A*) u ^ ^ 


T'- 
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Baldi et airi2O10l) ). together with other possi ble options for 
transfe r between DM and DE (see, e.g., lUopez Honorez et; 
( I2OI0I) ). The coupling 


C = b/m-p = ( 167 r/ 3 )^^^/ 3 /mp 


(3) 


sets the intensity of the energy flow from DM to DE. The op¬ 
tion of scale dependent C is also discussed in the literature, 
namely in connection with specific models (see, e.g., R. Mainini 
& S.A. Bonometto 2004). Let the background metric then read 

ds^ = a^{T){dT^ — dX^) , (4) 

(r : conformal time, dX : the spatial element); eqs. © then yields 

$ + 2--i> = -flV' -h Capc , Pc + 3-pc = -Cpc^ , (5) 
a a 

pc being the DM density, while V ($) is a self-interaction potential, 
as is required within quintessential DE models. 

The rational is to allow for an energy flow from (cold)-DM to 
DE. In this way, the field density could keep a steady fraction (some 
permits) of CDM density during the whole cosmic expansion, in 
spite of >1> being essentially kinetic above a suitable redshift z±. At 
z± the DE state parameter therefore shifts from ~ — 1 to ~ 4-1 
and such shift is found to be a generic feature, indipendently of the 
shape of V (4?). By keeping a significant field density at high z, this 
option eases one of the coincidence problems of ACDM. 


2.1 An early "L-CDM coupling 

Within a field theory context, a possible assumption is that CDM is 
a non-relativistic Dirac spinor field ip, interacting with $ through 
a Yukawa-like lagrangian 

Cm = -pf{^/m)ipip ; (6) 

here 2 mass scales, m = mpjb and p = grrip are introduced for 
dimensional reasons, nip being the Planck mass. In particular, b 
coincides with the coupling parameter in eq. (0 (see below). 

By assuming a kinetic part of the scalar field lagrangian Ck ~ 
dp$ dp^ , its equation of motion reads 

$ + 2--i-=-aV'-aVc4 : (7) 

a f 

once we work out that 

pc =-Cm = pf{b^/mp)ipip , (8) 

in the absence of a significant kinetic term for the spinor quanta. 
Notice that the number density operator for the spinor field n oc 
ipip, so that 

pc oc f{^/m)a~^ (9) 

according to the findings of lPas et"^ (l20Q6h . 

Eq. 0 is consistent with the former eq. 0 if/'// = -b/nip 
so that, if the function inserted in the Lagrangian 0 has the form 

/ = exp{—b^/nip) ( 10 ) 
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we re-obtain the coupled-DE theories of eqs. Also Pc accord¬ 
ing to eq. then exhibits a sort of exponential scaling, unless we 
assume that the argument of the exponential 


= ln(r/Tr) , 


( 11 ) 


Tr being a (unconstrained) reference time; then f = Trjr and 



( 12 ) 


Let us consider the radiative era, >1> bearing essentially kinetic en¬ 
ergy; then the r.h.s. of eq. 0 reads a^pcb/mp, while the whole 
equation means that pc = (rnplar')^ (thence pc oc a”'* !) or, 
equivalently. 


Stt 1 
3 m: 


2 (2/3 Pc) . 


(13) 


By comparing this equation with the Friedmann eq. , we deduce a 
constant early density parameter for CDM 


Oc 




(14) 


In turn, <1? being purely kinetic, owing to eq. GU, its energy density 


_ _ /mp\^ 1 _ 

“ 2a2 “ I, b ) 2a2r2 ^ 2 


(15) 


coinciding with its pressure. The state parameter of the future 
quintessential field is then w = -|- 1 , while its early density pa¬ 
rameter 


— 


Qc 1 
2 “ 4/32 ■ 


(16) 


In connection with the scale dependence found for pc, let us outline 
that, if we read Cm as a mass term of the Dirac spinor, its mass 
scales ex a~^ (as though being redshifted). As a matter of fact, 
both $ and CDM energy densities scale as a~'^. In the absence of 
coupling they would scale as and a~^, respectively. The flow 
of energy from CDM to $ is tuned to yield a faster (slower) CDM 
(4> field) dilution. 

The key issue, however, shown in BSLV12, is that these con¬ 
ditions are an attractor: even significantly perturbing “initial con¬ 
ditions”, the cosmic evolution rapidly settles on the behavior here 
described. Accordingly, such kind of evolution could last since a 
very early epoch, even since the late inflationary stages. 

We might tentatively assume that cosmic reheating was due to 
the very 4? field decaying into tp field quanta, a process stopping 
as soon as the attractor solution is attained. Interactions of t/) with 
other fields would then allow the 4? field energy to reheat other 
cosmic components, including both 7 ’s & v’s, as well as any other 
component then belonging to the thermal soup. 

The ensuing self-similar expansion, that can be defined 
conformally invariant (in the sense discussed, e.g., by Parker 
1969,1971), can only be broken if a component scaling differently 
from a“'* achieves a significant density. If the attractor solution is 
followed since inflation, 4> has had just a logarithmic time depen¬ 
dence since then. Its today’s value, therefore, is just ~ 60 times its 
value at the end of the reheating stages. 


3 EXIT FROM THE STATIONARY REGIME 

In any reasonable cosmological model, baryon density would even¬ 
tually break the conformally invariant expansion at a redshift 


Zb > 500. Adding just a baryon component to radiative compo¬ 
nents (coupled and uncoupled) is however insufficient to meet ob¬ 
servations. Among viable possibilities, BM14 outlined the option 
of including a Warm DM (WDM) component of thermal origin, 
with a temperature parameter Tw. Its derelativization occurs when 
Tw rriw (mass of WDM quanta) at a redshift Zw. At 2 ; S> Zw, 
being Pw — Pw/3, it is pw oc and WDM is one of the compo¬ 
nents of the thermal soup; then, at 2 2 w, being Pw ^ Pw and 

therefore pw oc Tw, WDM shall overcome the radiative component 
density so that the early stationary regime reaches an end. 

This also sets an end to the attractor solution of the coupled 
components. During their successive evolution, as in any model al¬ 
lowing for CDM-DE coupling with small or large /3, the state pa¬ 
rameter of the field component must turn from -l-l to ~ —1. This 
must occur about a suitable redshift z± so to allow a fair amount 
of today’s DE. In any approach based on a self interaction potential 
I/(4>), the fair 2 ± value is obtained by tuning some parameter(s) 
inside V (4?) expression itself. SCDEW models behave similarly. It 
is just convenient to make use of the first order field equation 

i l + u)„ 2 

4>i -b u)-4>i = —-—Cpca , (17) 

a 2 

instead of eq. 0; here 4?i = <!> and 2w = 1 -b 3iu — dlog(l -b 
w)/d log a. This equation, shown by BSLV12, requires w{a) to be 
assigned, instead of a potential 1^(4?). In this way, we select z± 
directly, without arguing about untestable V (T*) expressions. More 
explicitly, let us assume that 

w = ^^ with A = f — ^ (18) 

1 -bA \a±J 

with a± — (1 -b 2 ±)~^. In BM14 we also discussed the (mini¬ 
mal) dependence of results on the parameter e (= 2.9 here), whose 
arbitrariness replaces the choice of V (T*) expression. 

An option, not mentioned by BSLV12 and BM14, is that /3 
gradually fades, after the rise of p„ broke conformal invariance. 
No wanted feature of this class of cosmologies really depends on a 
late coupling and we shall see that, if /3 fading occurs late enough, 
just minor quantitative changes arise. 

Let us then consider the option that /3 oc exp(—a/a^g) with 
a-dg = Dx Ow. If we expect /3 fading to be a sort of consequence of 
the loss of conformal invariance, if should be D S> 1. Accordingly, 
different options will be labelled by the value of the delay (del.) 

d = logio39 (19) 

which is an extra parameter we introduce here. In Figure[T]we show 
the scale dependence of cosmic components in 4 different cases, 
passing from the option of ever lasting coupling, down to the case 
of delay d = 1. In this Figure and anywhere in paper I, we select 
the following parameter values: 

Do# Do 6 ha Tcmb ria 

0.7 0.045 0.685 2.726 3.04 0.968 

Here Do#. Dot, ha, Tcmb, Nv, rts are the present DE and baryon 
density parameters, the present Hubble parameter in units of 100 
(km/s)/Mpc, the CMB temperature, the number of (almost) mass¬ 
less i/’s, the primeval scalar fluctuation index, respectively. In the 
list, Dow and Doc are not included. Neglecting 7 ’s and u’s, their 
sum is 1 — Do# — Dob, while Do,c is fixed by the selection of the 
coupling constant /3 and Dow covers the residual gap. Notice that, at 
2 = 0, Do,c ■C 1/2/3^ being, typically, ~ 10~^Do(, (see Figure0. 

The option d = 4, that will be selected in Paper II, is char¬ 
acterized by adg — a±; as shown in Figure [T] for d = 4 density 
evolusion is just marginally affected by the fading of coupling. 
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Figure 1. Scale dependence of the densities of cosmic components after 
the break of conformal invariance. All models in this plot have /? = 10, 
mw = 90 eV. Delay (del) values from 1 to oo (ever lasting coupling) are 
considered. The thick curves labelled “cou” concern an ever lasting cou¬ 
pling. The behavior of densities does not suffer major changes, in respect to 
this option, unless d < 2. 


An earlier decoupling, instead, bears a number of conse¬ 
quences. Among them, a rise of today’s CDM density parameter 
flc,o, because the energy leaking from CDM to DE has a stop. If 
keeping d ^ 2, this rise is limited and, in this range of d values, the 
present CDM density parameter flc,o — 10“^ x fib,o' the present 
CDM contribution to the cosmic budget is 1/100 of baryons. 

A related effect is the formation of a dip in p<s>a* evolution. 
As a matter of fact, until $ is purely kinetic, its density keeps 
p<f = mpl(2ba^T^), thanks to the energy inflow from CDM. In 
the radiative era, when r oc a, is then constant. The transition 
to matter dominance should strengthen decline; however, when 
this occurs, coupled CDM density starts to increase, so providing a 
stronger <I> feeding. If the feeding stops before the kinetic-potential 
transition, occurring at , a progressive dilution of p* is unavoid¬ 
able. However, this produces a significant dip only if d < 2. Owing 
to these reasons, in the rest of this Paper I we shall never consider 
options d < 2. We shall also verify that models with any d > 2 
exhibit just minor differences. 

All through this discussion, we deliberately refrained from in¬ 
troducing any detailed physics causing a P fading. Possible mecha¬ 
nisms are discussed in Appendix A and further options surely exist, 
but the introduction of a further parameter can hardly be avoided. 
The treatment given in this work keeps on the phenomenological 
side. Just showing that /3 fading can be self consistently introduced, 
so allowing a plain treatment of nonlinear stages. 

Early nonlinear evolution will can be expected to induce a 
hierarchical rippled CDM distribution, made of “virialized” struc¬ 
tures, on growing scales. 

When and if /3 fades, CDM turns into a component of fast 
heavy particles. We expect them to cause no substantial effect on 
the linear evolution of other components, also because of the small¬ 
ness of the CDM density parameter; on this ground we decided to 
not consider these effects in our linear treatment. 

Before concluding this discussion on background evolution let 
us outline that the relation between expansion factor a (or redshift) 
and ordinary time t are almost identical in SCDEW and ACDM 


models; discrepancies never exceed 0.01% for any values of cou¬ 
pling /3 and delay d. 


4 LINEAR FLUCTUATION EVOLUTION 


In BM14 fluctuation evolution is studied in detail. Out of hori¬ 
zon initial conditions are determined and the system of differential 
equations, holding after the entry into the horizon, is numerically 
solved. Here, let us only remind some peculiar aspects, that this 
problem does not share with similar cosmological models. In par¬ 
ticular, field fluctuations are conveniently described by a variable 
(/5, related to the quintessential field /> as follows: 


(/ = $ H—, (20) 

mp 

$ being the background field. At the first perturbative order p> ful¬ 
fills the equation 


p + 2—/) + $/i + k\p + a V ($)‘p = 2/3 He ( — ) 5c • (21) 


Derivatives are taken in respect to the conformal time, as the metric 
reads 


ds^ = a^(r)[dr^ — (Sij + hij)dxidxj)] . (22) 


Moreover, Sc are the fluctuations in the coupled CDM component 
and h, the trace of the 3-tensor hij, describes gravity. Finally, k 
yields the mass scale M = (47r/3)p (27r/fc)^ of the fluctuation. 

Most terms in eq. (ED have therefore a transparent meaning, 
apart of the term U"((l?), where we apparently refer to a i/ self¬ 
interaction potential. Aiming to skip any reference to such quantity, 
BM14 show that one can use the relation 


2U" 


+ 


A f d e 

1+A\a1+A 

d $ d / a 
dr ya® 


$ 

eo + 2C^\ 

\ $ J 


-b 


(23) 


with 66 = e — 6. Here, A and e are the quantities defining the w{a) 
behavior in eq. ([HJ- All variables in this expression are known at 
each step of a linear evolutionary algorithm. 

It is then significant that, within this formulation, also in the 
fluctuation equations only 4? and 4? appear, while 4? values never 
matter. BM14 made use of these relations to modify the public 
algorithm CMBFAST and here we shall report results obtained with 
this code. 


4.1 Early growth of coupled CDM fluctuations 

Figure ID and [D then show a few examples of linear evolution, for 
fluctuations in the different cosmic components, from before their 
entry in the horizon, down to z = 0. The background cosmic pa¬ 
rameters are the same of Figure [D The mass of the warm com¬ 
ponent is anywhere 90 eV, its range being fixed by the arguments 
discussed below, in the next subsection. The k values selected for 
Figure ID where the delay parameter is d = 2, correspond to the 
mass scales of 1.8 x lO^'^/i”^ M© (or 8 ft“^Mpc , top panel), 
1.4 X Mq (central panel), and 8.7 x M© (lower 

panel). Figure |D is for 2 intermediate scales and d — 3. For the 
sake of comparison, the evolution of baryon and CDM fluctuations 
in a ACDM model with the same cosmological parameters is also 
overplotted. 
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Log a 



Log a 


Figure 2. Evolution of density fluctuations in the model indicated in the 
frame for 3 k values; from top to bottom they correspond to masses 
1.8 X 10^^- 1.35 X 10^^- 8.7 x respectively. The first 

k value then corresponds to a comoving length scale of 8 /i~^Mpc, at the 
boundaries between linear and non-linear behaviors; the second and third 
k values lay at the top and bottom limits of the galactic mass scales. Dotted 
curves show the evolution of CDM and baryons, for the same scales, in a 
ACDM model. 



Log a 



Log a 

Figure 3. As previous Figure, for d = 3, and 2 intermediate scales. 


A first result is shown in the top Figure the low -2 growth 
of fluctuations, namely on (scales close to) the linear regime, is 
the same in ACDM and SCDEW cosmologies. Numerical outputs 
confirm that discrepancies between ACDM and SCDEW growth 
factors, for k < 0.5, are below the expected precision of the algo¬ 
rithms used. It is so in spite of the (slightly) different final ampli¬ 
tudes. In the non-linear regime, discrepancies are appreciable, al¬ 
though still small. Aside of that, it is worth outlining that SCDEW 
models are also indistinguishable from ACDM, as far as CMB 
spectra are concerned. Discrepancies between SCDEW and ACDM 
CMB spectra, for /? > 5 are again below the expected precision of 
the algorithms used. 

Another evident feature is the fast growth of coupled CDM 
fluctuations, not discontinued at the time when the horizon attains 
the fluctuation size. In ACDM, uncoupled CDM fluctuation grow 
much more slowly, namely until radiation exceeds CDM density, as 
radiation fluctuations (still coupled to baryons) have stopped their 
growth, being in the sonic regime. In SCDEW models, this extra 
growth is expected, in spite of the small density pc ~ pt/2/I^ 
(pt : total density) of the coupled CDM component: in the non- 
relativistic regime, coupling causes an increa se of the effectiv e 
CDM-CDM gravity by a factor 1 + (4/3)/3^ ( lAmendolaHl999h : 
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for S> 3/4, this factor erases the division by 20^ so that, as 
far as self-gravity is concerned, it is as thought CDM had a density 
close to the total density pt- In addition, CDM particle velocities 
are enhanced by a continuous push due to their progressive mass 
decrease (the so-called antidrag term). All these effects are taken 
into account by the linear program and cause an increase of CDM 
density fluctuations, after their entry in the horizon and in the radi¬ 
ation dominated epoch, approximately oc . 


4.2 WDM fluctuation regeneration 

Large amplitude CDM waves cause a restart of fluctuations in 
WDM, as soon as it becomes non-relativistic. As is also visible 
in Figures 12 and 0 WDM derelativizes at Zw.der ~ 10®. Derela- 
tivization is linked to the value taken by the T^/T-y ratio, between 
WDM and photon temperatures. If WDM particles decoupled from 
the thermal soup when its density was /30)g„^decT'^, today 
(T„/r.y)® ~ (2/grw,dec)(l + ‘iN^/22) : number of (almost) 

massless i/’s in the cosmic background; gv,,dec' number of relativis¬ 
tic degrees of freedom at the time of WDM decoupling), while 


4.3 WDM particle mass selection 

WDM fluctuation regeneration in SCDEW models contrasts with 
undisputed free streaming effects in standard AWDM models, 
causing a large-fc cutoff to the spectral function 

A"(fc) = ^k^Pik) , (26) 

(P{k) = (|5(fc)|^) : transfered spectrum), because of the erasing of 
any fluctuation entering the horizon before WDM has derelativized. 
With mw 90 eV, the decline of starts at fc ~ 0.5 LMpc”^, 
being already damped by a factor ~ 10® at fc ~ 1 fcMpc“®. Ac¬ 
cordingly, viable AWDM models currently refer to WDM masses 
rriw — 2-3 keV, allowing fluctuations to survive up to fc ~ 20- 
30 /iMpc"^ 

Unfortunately, the need of such a large m„ partially invali¬ 
dates the choice of warm instead of cold DM. The residual mo¬ 
tions of low-mass particles, e.g., are then insufficient to reduce the 
number of expected MW satellites or to prevent them settling on a 
NFW (Navarro, Frank & White 1997) profile. Recent simulations 
c onfirm a relation betw een m„ and the size of a core. According 
to lMaccio et al.l ( l2012h . in the dwarf galaxy mass range, the size of 
the core 




45C(3) mw 

47r-‘ \TjJ lO^Tj ' 


(24) 


For Ni, = 3 we have then 


^ 0 . 115 ’^"'^^^ (25) 

Qw,dec 

so that T^jT.y ~ 0.32 if — 0.12 and mw — 90 eV. 

Accordingly, T„ ~ mw at a redshift ~ 10®-10® , in agreement 
with what is shown in the plots, these effects being all included in 
modified CMBFAST. In particular, our algorithm takes into account 
the contributions of coupled CDM and >1> field to the thermal soup. 

At the time of WDM derelativization, CDM fluctuations Sc 
exceed then any other fluctuations, by a scale dependent factor 
Fb ~ 10^-10® for galactic scales. There is no boost to the CDM- 
WDM gravity due to coupling, but the CDM density excess Spc = 
pcSc — {pt/2,0^)5c is boosted by the factor Fb, compensating the 
division by 2/3^ ~ 10^ due to the small CDM density. As the fac¬ 
tor Fb depends on the time elapsed since the entry in the horizon, 
the restart of WDM fluctuations should be more and more effective 
towards greater fc values. In turn, WDM particle escape velocity 
from smaller size fluctuations is smaller, so that the effect is only 
partially visible in the final spectra. 

On the contrary, this very effect is unchallenged in baryons, 
whose mean velocities are negligible. Also the (later) restart of 
their fluctuations is mostly due to CDM. In the absence of residual 
particle velocities, baryon fluctuations growth even exceeds WDM 
fluctuations, both at high z and at 2 : = 0, namely at high fc. 

Since ~ ^w.der. WDM density starts to dilute just oc a“®, 
and primeval conformal invariance is broken. Also baryon density 
growth violates such invariance, of course, but baryon density over¬ 
comes radiative components at a redshift < 2w,tier- In principle, it 
can make sense that a break of conformal invariance and /3 fading 
are related events: when WDM has turned non relativistic, the /3 
coupling appears like a residual pleonasm. Leaving apart detailed 
options to model a relation between these effects, we just quantify 
it with the parameter defined in eq. 01. For instance, for d = 2 
(3, 4), P fades when 2 ~ 5000 (500, 50). Values of d 3> 6 return 
an ever lasting /(-coupling. 


R 


core 


/ 1.0 \ / lOOeV y' 

\/i“ikpc/ \ mw / 


(27) 


If mw ~ 2keV, e.g., it is rcore ~ 5/i“^pc. On the contrary, a 
fair core size, fitting observations, apparently requires mw ~ SO¬ 
HO eV; a mass scale yielding fair cores but no galaxies. It is then 
significant that the linear theory of strongly coupled-DE models 
exhibits a full restart of fluctuations. 

Before concluding that a solution to low-scale ACDM prob¬ 
lems is found, however, we shall first consider in detail the full 
shape of linear spectra and use them to perform ad-hoc numerical 
simulations. The former aim is fulfilled here below, the latter one is 
the target of Paper IT 

Meanwhile, let us outline that there are a few other elements 
in favor of choosing mw — 90 eV, indipendently of the expected 
halo and satellite predictions. The very Eigure[T]shows that the con¬ 
tributions of WDM, CDM and 4? to the primeval thermal soup are 
all within half order of magnitude. Even though this is due to the 
choice of P and mw, it would be reckless to conclude that it fixes 
their ranges, but it would also be hard to believe all that to be fully 
casual. 

Eigure|4]then shows a strictly related point: The primeval den¬ 
sity parameters Hi plotted there are for CDM, <!> and WDM. The 
WDM density parameter exhibits just quite a mild P dependence. In 
fact, only for very low P values, its level is significantly eaten by a 
non-negligible ’SjAp'^ contribute. On the contrary, the WDM level 
would siginificantly depend on mw. Should this mass be greater by 
Am, WDM derelativize earlier by Aa ~ Am/To and the WDM 
line should be lower. In order to recover the nearly-coincidence be¬ 
tween primeval Hi, fairly greater P values should then be selected, 
as the dependence on P is quadratic. The opposite would occur for 
lower mw values, with the risk to approach the limits of the physi¬ 
cal range for /(’s. 

This plot should be however taken together with a third ar¬ 
gument, which can be made only after giving details on spectra 
(see Eigure[2 below): equal normalizations for CMB angular spec¬ 
tra and as, in ACDM and SCDEW models, should be required; 
this fixes a narrow P range, once mw = 90 eV is selected. In 
turn. Figure |4] returns mw = 90 eV for that P range. Altogether, 
the two constraints are close to a system of 2 algebric equations 
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Figure 4. Primeval density parameter of WDM and coupled CDM-4? com¬ 
ponents in the radiation dominated era, for rriw = 90 eV. The upper (lower) 
dashed line yields CDM (<I>) density. The total density of the coupled com¬ 
ponents is the “parallel” solid curve. The slightly rising solid line is WDM 
density parameter. 


with 2 unknown, although yielding softer constraints. /3 ~ 10 and 
rriw = 90 eV are at the middle of the narrow allowed bands. If 
N-body simulations will confirm that this selection favours obser¬ 
vational profiles, etc., we may say that this follows from the optimal 
parameter choice within the model. 


4.4 CDM non-linearities 

Recovering (WDM or) baryon fluctuations thanks to coupled CDM 
action, on those scales where they had been erased, requires CDM 
fluctuations 100-1000 times wider than uncoupled CDM fluctua¬ 
tions in ACDM models. In fact, coupled CDM has a much smaller 
density parameter than usual CDM, when its action is needed. 
Luckily enough, such wide amplitude is an intrinsical feature of 
SCDEW models. In turn, this feature can cause the birth of early 
CDM nonlinearities. These non-linearities are anyway suppressed 
by the fact that CDM contributes to the cosmic budget as little as 
~ 1 per cent of the very baryonic component. 

In SCDEW models early non linear CDM structures how¬ 
ever form on small scales. In the radiation dominated period, after 
entering the horizon, coupled CDM fluctuations grow oc~ . 
If Jc ~ 10~® at the horizon redshift Zh, it shall be (5c ~ 1 at 
Zn.i. ~ X Accordingly, fluctuations on scales entering 

the horizon at a redshift 10 ® times greater than equality, involv¬ 
ing a mass ^ iCh~^M q are quite likely to produce CDM non 
linear structures. Higher mass structures can also form up to a scale 
depending on the delay d. 

The most immediate problem related to CDM non-linearities 
concerns the reliability of linear algorithms predicting fluctuation 
evolution. In fact, the first effect due to non linearity onset, is an 
acceleration of the growth rate. An estimate of the effect can be 
desumed from the theory of spherical growth in pure CDM mod¬ 
els, telling us that a density contrast ~ 5-6 is attained when the 
linear 5^ — 1.0-1.1 . Non-linearity triggers other effects, as mode 
mixing, a loss of Gaussianety and high velocity fields. Eorming 
structures, in particular, shall have a shrinking radius and a global 


motion, so that we can expect a substantial loss of coherence be¬ 
tween the distribution of CDM and other cosmic components, with 
top spectral contribution eventually affecting wavelengths smaller 
than the original fluctuations 5c- Accordingly, no excessive effects 
are expected on scales never overcoming a mild non-linearity. 

At an advanced stage, evolved CDM non-linearities shall then 
be non-linear structures embedded in an almost unperturbed con¬ 
tinuum, taking also into account that, altogether, the CDM density 
parameter has reduced to (!?( 10 “®- 10 ~'*), so that the perturbing 
bodies shall be relatively rare. 

Let us also outline that, also after f) fading, the linear algo¬ 
rithm will continue to treat CDM and baryons separately. As a con¬ 
sequence, the CDM fluctuation 5c can still overcome unity while 
the baryon fluctuation 5t, is still <C 1. However, in these conditions, 
CDM and baryon obey the same equations of motion and the phys¬ 
ical variable is 


^ _ Hc(5c + ^b5h 


(28) 


Unless this is non-linear, physical non-linearities do not exist. Also 
if, slightly before /3 fading, 5c approached unity, but without badly 
loosing coherence, the physical post-fading variable will be Scb', 
however, even if coherence is partially lost, that part of 5c still co¬ 
herent with 5b will enter a unified growth regime. Accordingly, the 
effect of marginal non linearities can be thought to be just a slight 
boost of the CDM-l-baryon fluctuation amplitudes. A quantitative 
estimate of such boost is not immediately feasible, although one 
must however recall that flc -C tdb, so that the very uncertainty 
shall not be too large. 

When considering the spectral function in the next Section, 
which are a product of linear algorithms, one must therefore dis¬ 
criminate among different possible meanings of apparently high 
A^{k) for the CDM component at high k. 

There is however a clear conclusion for the above discussion: 
SCDEW models predict the formation of early non-linear CDM 
structures up to a mass ~ 10^-10®With a present CDM 
density ~ lO®/i~^M 0 /(/i“^Mpc)®, we expect their average co¬ 
moving distance to be ~ 1 ft“'^Mpc, a value comparable with the 
average intergalactic distance. Smaller residual CDM structures 
could also exist, if not embedded in greater structures in a hierarchi¬ 
cal formation process. No further speculation is however possible, 
to better define the specific features these “objects” could exhibit, 
apart of the rather obvious but probably simplistic statement that 
they might occupy the nucleus of existing galaxies. 


5 FLUCTUATION SPECTRA 

The modified version of CMBFAST enables us to evaluate the trans¬ 
fer functions and the spectral functions A'^{k,z) (see eq. I26l l. 
Eigures show them at z — 0, for specific models, with sepa¬ 
rate curves for different components. Results are compared with 
ACDM. 

Eor log fe — 1 the spectra of all cosmic components almost 
overlap. At greater k values, baryons and WDM exhibit a gap in 
respect to ACDM. In the absence of strongly coupled CDM, such 
hint to decline would turn into a fast decrease at slightly higher 
k’s, a characteristic of standard AWDM cosmologies. Here, on the 
contrary, thanks to the action of coupled CDM, we see a fast recov¬ 
ery. Baryon and WDM spectra therefore re-approach the ACDM 
spectrum and, in the case d = 4, overcome it. In turn, the (cou¬ 
pled) CDM spectra are constantly above baryon and WDM; for 
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Log(k/Mpc"‘h) Log(k/Mpc'‘h) 


Figure 5. Spectral functions obtained from the linear theory at 2 : = 0 for models with d = 2 and 4. We compare them with ACDM spectra, finding that the 
revival of high-A: fluctuations is highly efficient. 



Log(k/Mpc"'h) 


Figure 6. High-fc difference between baryon and WDM transfer func¬ 
tions T(k), shown to be independent from d (we overlap curves for d = 
1.8, 2.4, 3, 3.5, 4). The ratios T 1 /T 2 are between the transfer functions 
for baryon or WDM and the total transfer function, as indicated in the Fig¬ 
ure. The model is m„ = 90 eV and /3 = 10. As expected, the split starts 
at greater (smaller) fc if a greater (smaller) m,„ is taken. 

d = 2, however, the CDM spectrum does first decline, as baryons 
and WDM do, then gradually rises above them; on the contrary, 
for d = 4, the CDM spectrum keeps constantly above ACDM and 
there is just a hint of slower increase where baryons and WDM 
begin their gap. 

The most significant point, however, is the comparison with 
ACDM. Figures show a greater spectral discrepancy for lower d 
values. If we however consider P ^ 10 values (not shown in the 
Figures), we have a similar evolution of CDM fluctuations, for the 
P compensation discussed in Section 4.1, but equally large CDM 
fluctuations are less efficient to recreate WDM and baryon fluctua- 



- 2-10 1 2 

Log R/h"‘Mpc 


Figure 7. Linear mass variance on the scale R for the model in the frame 
(an exponential window is used). 


tions, because of the smaller CDM mass attracting them. Accord¬ 
ingly, the discrepancy with ACDM is even greater, in spite of CDM 
having spectra similar to the ones shown here. As a matter of fact, 
with P — 10, the discrepancy of baryon and WDM spectra with 
ACDM is not so great, never exceeding 1 order of magnitude. 

There is another point, already visible in the spectral func¬ 
tions in Figure but conveniently stressed in Figure [b] at high k 
the baryon spectrum starts to exceed the WDM spectrum. It is not 
a negligible effect. At A: ~ 30feMpc“^, the baryon spectrum ex¬ 
ceeds average by ~ 60 %, while tha WDM spectrum exhibits a de¬ 
ficiency, still in respect to average, by ~ 25 %. In Figurecurves 
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referring to models with d = 2, 3, 4, are shown to overlap. The 
Figure is done for rriw = 90 eV; as expected, a significant increase 
of TTiw reduces the gap. 

It may be significant to see the shift between SCDEW and 
ACDM, when looking at the linear m.s. fluctuation 

dk A^{k) W^{kR) (29) 

(mass variance). There are different forms for the window func¬ 
tion W(kR)\ in this work all computations were coherently per¬ 
formed by using an exponential window. The (tr behavior is shown 
in Figure 13 for d = 4 where, within the models considered here, 
it is greatest. As recalled in the frame, this is a linear computa¬ 
tion, and the slight discrepancies appear on scales R < 1 /i“^Mpc 
(~ 3.5 X Mq), surely non-linear today. The expected 

overall effect, after the extensive non-linear fc-mode mixing, is ex¬ 
pected to be slightly greater. 

The expression J29t can be also used to study the /3 depen¬ 
dence of as (mass variance on the scale of 8 d~^Mpc) on /3 and 
rriw. In Figure [S] we keep rrim fixed and study the as dependence 
on P for models with d = 4. The main finding shown in this plot is 
that, in order to approach the as value of ACDM, we need ~ 10. 

Let us then outline the difference between total spectra of 
models differing for the d values. As shown in Figure[3 discrepan¬ 
cies increase with k. For logic ~ 1.5 (M ~ 2.8 x Itf h~^M q) 
the discrepancy is ~ 2, when we compare d = 4 and 2 ; it be¬ 
comes < 20 % when we compare d=3 and 4. The above feature 
however exhibits just a mild dependence (not shown in the Fig¬ 
ure): for lower (higher) /3 values, the discrepancy is slightly smaller 
(greater). 

These spectra were accurately considered to choose the 
SCDEW model most suitable to study non-linear effects. Coupled 
CDM non-linearities are safely milder with d = 2. (They would 
be even milder for d < 2 but, as shown in Figure 1, the d de¬ 
pendence of the model becomes stronger and, in particular, a sig¬ 
nificant gap in p<j behavior becomes unavoidable.) As discussed 
in previous Section, even in the d = 2 case we expect some im¬ 
pact of forming CDM non-linearities, probably slightly rising the 
spectrum, more and more as k increases. Owing then to the limited 
difference between spectra, shown in Figure|3 we decided to make 
use of the d = 4 spectrum, assuming that, probably, it best repre¬ 
sents the d = 2 case, once the upward shifts due to forming CDM 
non-linearities are included. 



5 10 15 20 



/S 


Figure 8. Lower plot: crs dependence on /3 for SC models best fitting CMB 
fluctuations (hence, normalized as shown in the top plot). In the SC model 
selected rriw = 90 eV. Dashed horizontal lines yield normalization and as 
for a ACDM model with the same parameters. 



Log(k/Mpc 'h) 


6 DISCUSSION 

A number of experiments are in progress, aiming to inject the 
nature of dark cosmic components, e.g. BOS^l HETDEjfl DEfl 
LSqSlfland Euclitfl Even assuming a linear expression w = wo + 
uii(l — o) for DE state parameter, however, errors on wo and 
can hardly be pushed bel ow 10 % and some 10 %, respectively (e.g. 
Iloachimi & Bridlj|2010l) . Better hopes exist to detect violations of 
the standard relation between scale factor a and growth factor Q 
time dependences. 

They could arise from violation of General Relativity if, e.g.. 


^ http://www.sdss3.org/surveys/boss.php 
^ http://hetdex.org/hetdex 
® http://www.darkenergysurvey.org 
^ http://www.lsst.org 
5 www.euclid-ec.org/ 


Figure 9. Dependence of the overall spectrum on the d parameter. At fc ~ 
32/iMpc“^ (M ~ 2.8 X 10^h~^M q) the ratio between d = 4 and 
d = 2 spectra is < 2 . 


the gravitational action R (Riemann scalar) is replaced by a suit¬ 
able function f{R)- In this case the dynamics of any component, 
baryons and DM, would be directly modified at large distances. A 
large deal of work has deepened t his option. N-body si mulations 
have been also recently performed dPuchwein et al.ll2013h . 

An alternative option is that only DM dynamics is modified, 
as occurs in “ordinary” coupled DE theories. Also this option has 
been widely debated, N-body sim ulations of these models have 
been performed since sev e ral yea rs i Maccio et al.l2004l : [Baldi et^ 
120101 : lui & Barrowll201 ll : [Baldill2012 ). finding that baryon distri¬ 
bution would also change, over very large scales, as an indirect 
consequence. Since a few years these models have been shown to 
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be consistent with data if /^ < 0.1 , limits being widened if cou¬ 
pling is considered together with a non-vanishing neutrino mass, 
even finding a ^ 2-(j constraint on the coup l ing 13 about ~ 
0.09 jLa Vacca et alll200^ : iKristiansen et^bOld : IPettorino et al.l 
l2012h . More recently, the coupling option was also shown to ease 
the ten sion betwee n Planck and Hubble telescope Hq estimates, al¬ 
lowing IXial 120 13h to turn limits into a > 3 -(t detection, yielding 
P = 0.078±0.022 and, consistently. Ho = 74.8±2.8. Also within 
this option, significantly affecting the dynamics of the whole DE, 
future data are expected to outline apparent a{t)-Q(t) discrepan¬ 
cies. 

Strongly Coupled cosmologies (SCDEW models), on the con¬ 
trary, predict no change in a{t) or Q{t) in respect to ACDM and 
therefore, a fortiori, the a{t)-Q{t) relation should be found do be 
consistent with ACDM predictions. The very equation of state of 
DE remains somehow arbitrary, although state parameters m < — 1 
would require suitable extentions of the approach described in this 
work. All above experiments are therefore expected to yield results 
seemingly consistent with ACDM cosmologies. 

This point has been accurately verified in this work. In com¬ 
parison with BSLV12 and BM14, two further points were outlined 
here: (i) By providing a Lagrangian approach to the coupling be¬ 
tween d? and tp fields (DE and CDM, respectively) and stressing 
the logarithmic growth of the former one, we outlined the possibil¬ 
ity that d> is both inflaton and DE. (ii) The risk of CDM fluctuations 
reaching an early non-linear dynamics was also outlined and cir¬ 
cumscribed. 

We shall further discuss the (i) point in another work. As far 
as the (ii) point is concerned, it is strictly linked to the minimal as¬ 
sumption of a constant interaction constant jd. We tentatively sug¬ 
gested to overcome such option by admitting fi to fade, once it 
accomplished its aim to allow us an (almost) conformally invari¬ 
ant cosmic expansion through cosmic ages. Possible mechanisms 
directly relating conformal invariance break to fading, however, 
were not discussed. A number of alternative options are briefly dis¬ 
cussed in Appendix A. 

Altogether, the picture considered remains fully viable and 
leads to a phenomenological picture quite close to ACDM. In its 
theoretical framework, however, most of the unpleasant fine tun¬ 
ings and coincidences of ACDM are significantly eased. 

SCDEW and ACDM cosmologies are however distinguish¬ 
able through a number of observables; First of all, SCDEW sug¬ 
gests a DM particle with mass ~ 100 eV. If this kind of warm-hot 
DM replaces the CDM of ACDM models, in the absence of strong 
coupling there would be almost no structure in the Universe. Even 
the option of mixed cold and hot-warm DM is far from fitting data 
on dwarf rotation curves or large galaxy satellites. Accordingly, dis¬ 
covering a sterile neutrino or a gravitino in the 100 eVmass range 
would mark a strong point in favour of SCDEW models. 

The other critical prediction of SCDEW cosmologies is the 
formation of early coupled-CDM structures: Coupled-CDM fluc¬ 
tuations Sc, to be able to revive baryon and WDM fluctuations, 
had to be quite large when WDM particles finally became non- 
relativistic. This is what theory predicts and is also the only way 
to make their density excess Spc = PcSc significant; in fact pc — 
\l2fP' or less. Therefore, Sc approaches unity earlier than other 
components. If, meanwhile, the CDM-ifi coupling fades, baryons 
and CDM share identical equations of motion and only their total 
fluctuation matters. Non-linearity then becomes just a formal and 
harmless feature. Should it be not so, bound CDM structures form. 
This however occurs over mass scales < 10^h~^M q (let us how¬ 


ever recall that, altogether, coupled CDM density is 1 % of the 
very baryon density, today). 

This risks to weaken the predictivity of linear codes. When 
growing non linear, CDM fluctuation gravitation becomes stronger 
than linear codes compute, speeding up the growth of other com¬ 
ponents. This effect is however weakened by pc becoming even 
smaller than ~ 1120^, as the Universe turns from radiation to 
matter dominated, and by the fact that other component fluctua¬ 
tions have recovered such a significant amplitude to yield an -at 
least- comparable push. Accordingly, we suggested the extra grav¬ 
itational push of CDM to bear effects equivalent to shifting the pa¬ 
rameter d from 2 to 4. Although based on qualitative arguments, 
this option is not unlikely to approach the real effect. Accordingly, 
N-body simulations in the accociated Paper II are based on the 
d = 4 option. 

All that implies a clear prediction: that early CDM primeval 
structures have formed. Although being extremely rare because of 
the low overall CDM density (~ 1/100 of baryon density), they 
could be the main existing structures at high redshifts {z ';$> 10) 
and might have an effect at later times on the formation of the first 
stars and/or cosmic reionization; we plan to address these issues in 
future work. 

In our opinion, however, the strongest argument in favour of 
SCDEW models is the apparent easing of fine tunings. In particular, 
a twofold DM is an option often considered on purely phenomeno¬ 
logical bases. However, the true DM, in SCDEW cosmologies, is 
the warm one; CDM is a sort of handyman component, fist allowing 
to keep on its tracker solution, then allowing inhomogeneities to 
revive on observational scales, eventually creating bound systems 
where it hides today. The real challenge of these models is finding 
out its hideouts. Apart of that, replacing A with a scalar field might 
seem not a fresh option. Also the independence from I.C., thanks 
to the presence of a field attractor, is not new. 

In respect to standard quintessence theories, however, 
SCDEW models are favoured due to a few strong points: (i) They 
are apparently independent from the choice of any specific poten¬ 
tial for the 4? field, (ii) The quintessential 4> field is shown to have 
necessarily been a significant cosmic component since long time, 
possibly since the end of inflation, (iii) Again, the attractor solution 
it fulfills is not linked to any potential choice, (iv) The possibility 
that DE and inflaton are the same field can also be pursued: in ef¬ 
fect, since inflation, 4? had just a logarithmic growth matching, e.g., 
the logaritmic evolution of constants in a Coleman-Weinberg-like 
potential, so able to provide potential energy both then and now. 

The idea that ACDM is such a successful! model because it 
mimics an underlyi ng more complex cosm ology is not completely 
new. For instance, lAmendola et alj ( l2008h investigated cosmolo¬ 
gies where DE arises because of the increase of neutrino mass, 
aiming to a scenario substantially indistinguishable from ACDM. 
Any such cosmology requires the introduction of extra parameters. 
In SCDEW models 2 extra parameters (apart of the delay d) are 
needed: the mass of the WDM particle and the coupling 0 How¬ 
ever, if we vary these parameters within a reasonable range, we ob¬ 
tain cosmologies which could also be mimiced by sorts of ACDM 
with a different choice of its basic parameters. Moreover, the mass 
m„ could soon become a non-free parameter, if a particle candi¬ 
date is found. In this case, the costraints on /?, deriving from fluc¬ 
tuation amplitude, CMB spectrum, density parameter choice, etc., 
as illustrated in this paper, are decisive. 
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APPENDIX A: PHYSICAL OPTIONS FOR /3 FADING 


In the CDM-<J> coupling C = b/mp one can straightforwardly 
introd uce a $ dependence, yielding Cf iMainini & Bonomettg 
( l2004h . e.g., took C = 6/$. In the context of SCDEW modes, ow¬ 
ing to the progressive increase, instead of taking /? oc {a/adg)°‘, 
we can assume /3 oc (<l?/$dg)“; an expression relating adg and a 
to <J>dg and a is then obtainable, so to achieve similar results. 


A smarter option, perhaps, still keeping to field theory, 
amounts to replacing the lagrangian ( 6 ) in Section 2 by 


Cm. = ■ 

This modifies the field and pc equations (5) into 


(Al) 


<i>-F 2 -<E- = -aV' +pc 


Ca^ 


1 -I- exp(C'4?)/I//i 


pc + 3-pc = -pc— - - , (A2) 

a 1 -F exp(C$)p/p 

where we keep to the expression of C = b/nip with constant b. 
Until 6 $ <C rup coupling is as usual. But, when 4? approaches the 
Planck scale, the effective coupling is cut off. 

A further option, not affecting field theory, can also be envis¬ 
aged, although its effectiveness depends on a number of still open 
problems. As a matter of fact one can wonder what happens to 
CDM nonlinearities after they entry in the horizon. The question 
is whether the nonlinear gravitational growth of CDM fluctuations 
ends up into a relativistic collapse or into virialized structures. In 
the former case, one should not forget that the action of on CDM 
can be described by replacing G with G*, adding the antidrag term, 
etc., only until a non relativistic approximation holds. On the con¬ 
trary, the dynamics of a relativistic collapse, in the presence of 4? 
interactions, is still an open problem. 

We can however formulate two conjectures: (i) The nonlinear 
gravitational growth ends up into virialization if the CDM matter 
density overcomes a suitable theshold pc. (ii) If relativistically col¬ 
lapsed CDM structures form, the <J> field is unable to interact with 
CDM there inside. 

Within this context, let us imagine to treat the growth of a 
spherical density fluctuation, starting from its entry in the horizon, 
when it involves fluctuations of amplitude Ac ~ A ~ 10~® (for 
CDM and the other cosmic components, respectively). Let R = ac 
be its radius, a being the scale factor. At the conformal initial time 
Ti, the radius is Ri = aid n with Ci = 0. Let then x = c/cp at 
any t > Ti,x shall fulfill the equation 


Ac'yfJc 


-F 2(1 - Dc)Aa 


(A3) 


almost exact in the period characterized by self-similar expansion. 
The equation is more complex later on, but the key features are the 
same. In particular, let us remind that Dc = 1/2/3^, 7 = 1-F4/3^/3. 

According to this equation, as expected, the CDM enhance¬ 
ment expands, enters a non-linear regime, and then starts to recon¬ 
tract, while the other components still undergo a linear evolution. 
Recontraction ends up into a complete gravitational collapse only 
if the Schwartzshild density 

p, ~ 3 X lO^®(g/cm®)(M/M0)-2 


is attained before arriving to the critical density pc. Here M is the 
mass involved in the collapse. Accordingly, at early times, CDM 
non linearities can be described as a rippled CDM distribution. 
There will however be a time after which CDM starts to form sorts 
of Black Holes (BH). 

In the late epochs, therefore, CDM has mostly turned into a 
suitable BH distribution, ceasing its interaction with $. 


This paper has been typeset from a TjiX/ ETjX file prepared by the 
author. 
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